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Abstract. Globular clusters as co Cen and NGC 2808 appear to have a population of very 
He-rich stars. From a theoretical point of view, one expects the presence of He-rich stars 
in all globular clusters showing an oxygen- sodium (0-Na) anticorrelation. In this paper, we 
briefly recall how fast rotating massive stars could be the main source of the material from 
which He-rich stars have formed. We speculate that the UV-upturn phenomenon observed 
in all elliptical galaxies might be due to He-rich stars. If this hypothesis is correct then fast 
rotating massive stars may have played in the early evolutionary phases of these systems a 
similar role as the one they played in the nascent globular clusters. 
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1. Chemical anomalies in globular 
clusters 

The finding of a double sequence injhe glob- 
ular cluster (jj Cen by I Anderson! (1 19971) and 
the further interpretation of the bluer sequence 
by a strong excess of helium constitutes a 
majo r enigma for stellar and galactic evolu- 
tion (lNorrisll2QQ4l) . The bluer sequence with a 
metallicity [Fe/H]= -1.2 or Z = 2 • 10"^ im- 
plies an He-content Y = 0.38 (0.35-0.45), i.e. 
an helium e nrichment AY = 0.14 (see also 
lNorrisll2QQ4 . In turn, this demands a relative 
helium to metal enrichment AY/AZ of the or- 
der of 70 (iPiotto et al.ll2005l) . This is enormous 
and more than one order of magnitude larger 
than the current value of AY/AZ = 4-5 
(iPagel et al.lll992h obtained from extragalactic 
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HII regions. A value of 4- 5 is consistent with 
the ch emical yields from supernovae (Maede3 
Il992h forming black holes ab ove about 20-25 
Mrr^. lMaeder & Mevnej (l2006h showed that the 
wind contributions of low Z rotating stars are 
able to produce the high AY/AZ observed in 
the blue Main Sequence in oj Cen. 

Recent observations show that He-rich 
stars in globular clusters might be a common 
feature. For instance, clusters like NGC 2808 , 
M13 and NGC 6441 (ICaloi & D' Antonal2007l) 
show a well-developed blue horizontal branch, 
and a strong slope upward from the red clump 
to the blue of the RR Lyrae regions. Both fea- 
tures could be explained if He-rich stars are 
present in these clusters. iKaviraj et al.l (l2007l) 
studied the UV and optical properties of 38 
massive globular clusters in M87 (elliptical 
galaxy in virgo cluster). A majority of these 
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clusters appear extremely bright in the far-UV. 
Canonical models (i.e. models without any su- 
per He-rich stars) would imply ages for these 
clusters about 3-5 Gyr larger than the age 
of the Universe according to the WMAP data! 
This difficulty is removed when a super-He- 
rich (AF/AZ >~ 90) stellar component is sup- 
posed to be present in these clusters. 

It has also long been known that globu- 
lar cluster stars present some striking anoma- 
lies in their content in other light elemental: 
while in all the Galactic globular clusters stud- 
ied so far one finds "normal" stars with de- 
tailed chemical composition similar to that 
of field stars of same metallicity (i.e., same 
[Fe/H]), one also observes numerous "anoma- 
lous" main sequence and red giant stars that are 
simultaneously deficient (to various degrees) 
in C, O, and Mg, and enr iched in N, Na, and 
Al (for recent revi ews see iGratton et al.]l2004l : 
ICharbonnell l2005l) . Additionally, the abun- 
dance of the fragile Li was found to be anticor- 
related with that of Na in turnoff^ stars in a cou- 
ple of globular clust ers dPasquini et al.l l2Q05l : 
iBonifacio e'taDl2007l) . 

It is clear now that these chemical peculiar- 
ities have been inherited at birth by the low- 
mass stars we observe today in globular clus- 
ters, and that their root cause is H-burning 
through the CNO-cycle and the NeNa- and 
MgAl-chains that occurred in an early genera- 
tion of more massive and faster evolving glob- 
ular cluster stars. In other words, compelling 
evidence leads us to believe that at least two 
generations of stars succeeded in all Galactic 
globular clusters during their infancy. The first 
one corresponds to the bulk of "normal" stars 
born with the pristine composition of the pro- 
tocluster gas; these objects are those with the 
highest O and Mg and the lowest He and Na 
and Al abundances also found in their field 
contemporaries. The second generation con- 
tains the stars born out of material polluted to 
various degrees by the ejecta of more massive 
stars, and which present lower O and Mg and 

^ On the contrary, the content in heavy ele- 
ments (i.e., Fe-group, of-elements) is fairly con- 
stant from star to star in any well- studied individual 
Galactic globular cluster (with the notable exception 
of CO Cen). 



higher He, Na and Al abundances than their 
first generation counterparts. 

Any model aiming at explaining the 
chemical properties of globular cluster stars 
should give an answer to the following ques- 
tions: (1) Which type of stars did produce 
the material enriched in H-burning prod- 
ucts? (2) What is the physical mechanism 
responsible for selecting only material bear- 
ing the signatures of H-processing? (3) Why 
does this process occur only i n glo bular 
clusters iR'Prantzos & Ch arbonnell (l2006h and 
Decressin et al. ( 2007b; a) propose answers to 
the above three questions. The main line of the 
reasoning is recalled below. 

2. Slow winds of fast rotating 
massive stars 

Fast rotating massive stars, with an initial 
metallicity Z = 0.0005, with typical time av- 
eraged velocity of 500 km s"^ on the main se- 
quence easily reach the critical velocit}0 at the 
beginning of their evolution and remain near 
the critical limit during the rest of the main 
sequen ce and part of the cent r al He-burning 
phase (iDecressin et al.l l2007bl: lEkstrom et al.l 
l2007h . As a consequence, during those phases 
they lose large amounts of material through a 
mechanical winfl which probably leads to the 
formation of a slow outflowing keplerian equa- 
torial disk. This is t he kind of process believ ed 
to occur in Be stars dPorter & Rivini us "2003"). 

The material ejected in the disk has two in- 
teresting characteristics: (1) due to rotational 
mixing that transports the products of the CNO 
cycle and of the Ne-Na and Mg-Al chains from 
the core to the stellar surface, it bears the sig- 
natures of H-burning and presents abundance 
patterns similar to the chemical anomalies ob- 
served in the second generation stars (see for 

^ Indeed up to now the peculiar chemical pat- 
terns observed in globular clusters, i.e., the 0-Na 
and Mg-Al anticorrelations, have not been found in 
field stars. 

^ The critical velocity is the velocity at the equa- 
tor such that the centrifugal acceleration balances 
the gravity. 

^ For example, our 60 M© model loses 24.3 M© in 
its slow wind. 
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more details iDecressin et al ] I2QQ71#I : (2) it 
is released into the circumstellar environment 
with a very low velocity and thus can easily be 
retained into the shallow potential well of the 
globular cluster. In the following we shall call 
it the slow wind. 

The disk- star configuration lasts during the 
whole main sequence for stars close to or at 
critical velocity; for the most massive stars that 
are at the QF-limit {i.e., when the surface lu- 
minosity is near the Eddington limit), it per- 
sists on the LBV phase. When the star evolves 
away from the critical limit due to heavy mass 
loss, the radiatively-driven fast wind takes over 
and the disk is not supplied by the star any- 
more. This happens during the central He- 
burning phase, before the He-burning products 
reach the stellar surface and contaminate the 
slow wind component (bottom panel of Fig. [T]). 
From that moment on the high-speed material 
ejected by the star and by the potential super- 
nova will escape the globular cluster. We pro- 
posed that this filtering mechanism, which re- 
tains in the potential well of the globular clus- 
ter only the slow stellar ejecta, is the phys- 
ical mechanism responsible for supplying the 
required H-processed material for forming the 
stars of second generation. 

The large star- to- star spread in light ele- 
ments observed today indicates that all the ma- 
terial ejected by the first generation massive 
stars did not have the time to be fully mixed 
before being recycled in the second generation. 
In the frame of our model, it imposes (see also 
Sects. 3 and 4.3), that star formation occurred 
in the vicinity of individual massive rotating 
star, from clumps made of slow wind material 
and of pristine interstellar gas. This star forma- 
tion may be triggered for instance by the ioni- 
sation front. 

The material ejected in the slow wind 
presents various degrees of enrichment in H- 
burning products. At the beginning of the disk- 
star phase, pristine material is ejected giving 
birth to stars with a composition similar to that 



^ We remind however that in order to reproduce 
the lowest Mg values observed we had to postulate 
an increase of the ^'^Mg(p,7)^^Al nuclear reaction 
rate with respect to the published values. 
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Fig. 1. Schematic view of the evolution of fast ro- 
tating massive stars. The colours reflect the chemical 
composition of the various stellar regions and of the 
disk, (top) During the main sequence, a slow out- 
flowing equatorial disk forms and dominates matter 
ejection with respect to radiative winds. Early on the 
main sequence, the [0/Na] value in the slow wind 
is the pristine one. (middle) As mixing proceeds, 
matter with [0/Na] typical of H-burning (i.e., ~ - 
2) is ejected in the slow wind (bottom) Due to heavy 
mass loss, the star moves away from critical velocity 
and does not supply its disk anymore; radiatively- 
driven fast wind takes over before the products of 
He-burning reach the stellar surface. Figure taken 
from .Decressin et al.. (.2007 a) . 



of the first generation. As time proceeds, the 
slow wind becomes more and more polluted 
in H-burning products (see Fig. 1 and ^ and 
forms stars that are more and more "anoma- 
lous". The slow wind phase lasts only a few 
million years {e.g., A. 5 Myr for a 60 M© star). 
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Fig. 2. Evolution of the surface abundances in mass 
fraction as a function of the remaining mass for a 60 
Mq stellar model at Z = 0.0005 and an initial ro- 
tational velocity on the ZAMS of 800 km s"^ The 
horizontal line labeled "Slow Winds" in the lower 
panel extends over the part of the evolution dur- 
ing which the rotating star releases slow equato- 
rial winds. Second generation stars are formed from 
clumps made of this material and of pristine inter- 
stellar medium. We can see how as a function of 
time the abundance of oxygen in the slow wind de- 
creases while that of sodium increases. 

3. Constraints on dilution factor, IMF 
slope and degree of evaporation 

How many massive fast rotators are needed 
in order to reproduce the observed popu- 
lation of anomalous stars in globular clus- 
ters? This question i s first addressed by 
iPrantzos & Charbonnell (l2006h for NG C 2808 
and in more details in lDecressin et ID (l2007ah 
on the basis of theoretical stellar yields for 
NGC 6752. Here we just recall the method 
used as well as the main results for NGC 6752. 

From rotating massive star models we can 
compute the amount of material ejected by 
slow winds with a given chemical composi- 



tion (see for details iDecressin et al ] l2QQ7bh . 
As already discussed in § 2, there is com- 
pelling evidence that the slow wind of an in- 
dividual star pollutes only on a small scale 
around its stellar progenitor, and that it is di- 
luted locally with pristine ISM. In order to con- 
strain quantitatively this local dilution process 
we use the Li vari a tion d etected in NGC 6752 
by IPasquini et al.l (l2005l) . It is important to 
recall that the massive star ejecta are totally 
"Li-free", as this fragile element is destroyed 
in these objects. We suppose that the mat- 
ter ejected by an individual massive star early 
on the main sequence encounters more pris- 
tine gas and is more diluted than the matter 
ejected later. This produces a Li-Na anticorre- 
lation. We chosed the dilution factor in order to 
re produce the Li-Na an ticorrelation observed 
by IPasquini et al.l (l2005l) . Doing so we obtain 
that the total mass of second generation stars 
amounts to about twice the total mass released 
under the form of slow winds by massive stars. 

Given a dilution factor, we can compute 
how many small mass stars can be formed from 
the slow wind of a given fast rotating massive 
star. With an Initial Mass Function (IMF), we 
can estimate the total number of second gener- 
ation low mass stars as well as the number of 
normal first generation low mass stars. Thus, 
the number ratio of anomalous to normal stars 
at the time of birth of the second generation 
stars can be predicted. During its evolution, a 
globular cluster may lose stars, for instance by 
tidal stripping. The ratio of anomalous to nor- 
mal star in that case will change as a function 
of time. Therefore, the number ratio of anoma- 
lous to normal stars in present day cluster de- 
pends on three parameters: the dilution factor, 
the slope of the IMF and the number of stars 
which have escaped the cluster (evaporation). 

As mentioned above, the dilution factor has 
been calibrated in order to rep roduce the Li- 
Na an ticorrelation observed by IPasquini et al.l 
(l2005l) . Its value only depends on the stellar 
models used. To constrain the IMF and the 
degree of evaporation we have only one ob- 
servable, the present day observed fraction of 
anomalous to normal stars in globular clusters 
(this v alue is equal to about 6 in NGC 6752, 
iCarretta et al.l l2QQ7l) . Thus we need to make 
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some additional hypotheses to constrain these 
two quantities. 

In scenario I we assume that the cluster un- 
derwent no evaporation of stars. In that case 
the present day observed fraction of anoma- 
lous to normal stars gives directly access to the 
fraction at the time when the second genera- 
tion was born. In this scenario the IMF can be 
constrained and the degree of evaporation is set 
equal to zero. 

In our Scenario II, we assume that the mas- 
sive polluters of first generation were born in 
the center of the cluster or have migrated very 
rapidly toward this region. Then the second 
generation stars are created only in the cen- 
tral region while the external part of the cluster 
hosts only first generation long-lived low-mass 
stars. This strong radial distribution is assumed 
to stay in place until the moment where su- 
pernovae sweep away the residual intra-cluster 
gas and thus strongly modify the cluster poten- 
tial well. Evaporation will then preferentially 
remove normal stars in the outer layers of the 
cluster. In the Scenario II, we impose a slope 
of the IMF equal to 1.35 and the degree of 
evaporation of normal stars is constrained (no 
anomalous stars are supposed to have escaped 
the cluster). 

To produce a high number of polluted stars 
(typically 5.7 times the number of normal stars 
as observed in NGC 6752) we need either (a) 
in Scenario I, a flat IMF with a slope for mas- 
sive stars around 0.55 or (b) in Scenario II, that 
about 96% of the unpolluted normal low-mass 
stars has been lost by the cluster. 

Using the two sets of free parameters 
obtained in Scenarios I and II, we then 
explored various consequences which are 
briefly recalled belo\y (see for more details 
iDecressin eraDl2QQ7ah . 

1) In both scenarios, we can produce an 
[0/Na] anticorrelation in the range -1 < 
[0/Na] < 0.62. Quite remarkably, this range 
corresponds to the one observed in NGC 6752 
if one considers only the set of stars for which 
both measurements of oxygen and sodium are 
available. Super-O-poor stars whose presence 
is indirectly inferred in some globular clusters 
need lower dilution coefficient or some addi- 
tional mixing occurring in the star itself. 




Fig. 3. Upper panel: Expected anticorrelation be- 
tween helium and [0/Na] for second generation 
stars (at birth!) born from material ejected through 
slow winds of different initial mass models. Lower 
panel: Same as upper panel for the correlation 
between ^^C/^^C and [0/Na]. Figure taken from 
IDecressin elaP (l2007ah 



2) In both scenarios, around 12% of the 
low-mass stars are expected to display a He 
content higher than 0.4 (in mass fraction). Note 
that in the present model. He-enrichments are 
obtained keeping the sum C-hN-hO equal to its 
initial value as requested by the observa tions 
(ICohen & Melendezll2QQ5l:ICarrettall2QQ6h . 

3) An anticorrelation is predicted between 
the helium abundance and the value of the ratio 
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[0/Na] at the surface of turnoff stars (see upper 
panel of Fig. [3]). 

4) More than 30% of stars are expected to 
have surface ^^C/^^C ratios between 10 and 30. 

5) The ratios ^^C/^^C and [O/Na] are corre- 
lated (see lower panel of Fig.[3]o. All the stars 
with [O/Na] ratios below about -0.6 are pre- 
dicted to show values of the ^^C/^^C ratios be- 
tween 15 and 20. 

6) Only a few percents of stars present at 
birth ^^O/^^O ratios of the order of a few hun- 
dreds. 

7) A correlation with a very large disper- 
sion is predicted between the ratios ^^O/^^O 
and [O/Na]. If only stars more massive than 
about 60 Mq would contribute to the building 
up of the chemical anomalies in globular clus- 
ters, then only values above 2000 would be ex- 
pected would be expected. 

8) A correlation between [C/N] and [O/Na] 
is expected, whose dispersion increases at 
lower [O/Na] values. 

9) In the case of Scenario I, the present 
day mass cluster represents about half of the 
mass which has been used to form stars. The 
rest consists in matter ejected by stars and lost 
by the cluster. The mass locked into remnants 
represents 19% of the mass of the present day 
cluster, the first and second generation stars re- 
spectively 12 and 69%. 

10) In the case of Scenario II, the present 
day mass cluster represents only about 9% of 
the mass which has been used to form stars. 
The rest consists in matter ejected by stars and 
of stars and stellar remnants lost by the clus- 
ter. The mass locked into remnants represents 
8% of the mass of the present day cluster, the 
first and second generation stars respectively 
16 and 76%. 

The above estimates are based on observed 
caracteristics of NGC 6752. Using other clus- 
ters for calibrating the free parameters would 
of course produce diff'erent results. 



^ These correlations are obtained for non-evolved 
stars. When the star evolve along the red giant 
branch, internal mixing may change somewhat the 
results. 



4. Discussion 

4.1. Why no anomalous stars in the field 
of the halo? 

Globular clusters contain a significant popu- 
lation of anomalous stars. In the field, at the 
moment, only normal stars are observed at 
the corresponding metallicity (no anticorrela- 
tion found). This indicates that something spe- 
cial occurs in clusters which does not occur in 
the field. What might it be? There are many 
possibilities: 1) stars in the field were born in 
much shallower clusters which were rapidly 
destroyed when the first supernovae appeared. 
This disruption process prevents a second gen- 
eration of stars to appear; 2) in small easily 
destroyed clusters, too few massive stars are 
present, thus the source of CNO-rich material 
is too weak for allowing many anomalous stars 
to be formed; 3) another possibility would be 
that the rotation rate of massive stars depends 
somewhat on the strength of the star formation 
event. This is of course quite speculative, but 
it is interesting to note that the rotation rate of 
B-type stars in clusters seems to be higher than 
that o f B-type stars in the field (iHuang & GiesI 
l2006h . 

Even if anomalous stars are formed exclu- 
sively in clusters, one can still wonder why 
some of them could not have escaped the clus- 
ter and passed the rest of their life in the field. 
It is reasonable to think that during their life- 
time the globular clusters lose part of their stel- 
lar content by "evaporation" processes like for 
instance tidal disruption. Thus, why no anoma- 
lous stars are lost in that way? As explained in 
our scenario II, this may arise due to the fact 
that anomalous stars are preferentially formed 
in the central parts of the cluster and are still in 
that region when the major evaporating events, 
aff'ecting mostly the stellar population lying in 
the outer layers of the cluster, occur. 

It might be also that anomalous stars are 
actually present in the field, but in such a small 
proportion that they have until now escaped de- 
tection. Let us estimate an upper limit to the 
fraction of anomalous stars in the field, /a. If 
no anomalous stars are detected among A^obs 
field stars, it means that ^A^obs « 1 or that 
/a « 1/A^obs- In case A^obs ~ 100, much 
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less than one star over hundred in the field is 
anomalous. 

From A^obs and the degree of evaporation 
(see below), we can constrain Fa, the frac- 
tion of stars lost by globular clusters which are 
anomalous. Let us recall that the number of 
stars in the halo is around 10^. The number of 
stars in present-day globular cluster amounts 
to ab out 1 % of the total number of stars in the 
halo (IWoltiedll975b . Let us say that the glob- 
ular clusters have lost K times the number of 
stars they contain today. This implies that the 
fraction of stars now in the field but having 
been born in globular clusters is / = /sTO.Ol. 
If ^ = 10 (i.e. if globular clusters lost about 
10 times the number of stars they contain to- 
day), then one field star over 10 would have 
been born in a globular cluster. Then the frac- 
tion of anomalous stars in the field, fA, can 
be written ^ = « 1/A^obs- This im- 

plies that Fa « With K = 10, and 

A^obs ~ 100, Fa « 0.10. Increasing A^obs will 
provide tighter constraint on Fa (provided K is 
known from other considerations). From A^obs, 
an IMF, as well as an observed averaged frac- 
tion of anomalous to normal stars in globular 
clusters some constrains on the degree of evap- 
oration of first and second generation stars can 
be obtained. 

4.2. Supernova rates in very young 
globular clusters 

In case the chemical anomalies are linked to 
a population of massive stars, what would be 
the consequences on the expected rate of super- 
novae in the early evolutionary phases of glob- 
ular clusters? With a slope of the IMF equal 
to 0.55 as in our scenario II, one would ex- 
pect a supernova rate of 1 to 2 per 1000 year 
per globular cluster (containing initially one 
million stars with masses distributed between 
0.01 and 120 M©). This rate would be main- 
tained for a few million years (about 8 mil- 
lions years). If the halo of a galaxy contains 
about 100 globular clusters, then, one would 
expect 1 to 2 supernovae per 10 years observ- 
ing the system composed by all globular clus- 
ters. Interestingly, part of these supernovae. 



those originating from initial masses superior 
to about 50 M© would likely give birth to col- 
lapsars and then to Gamma Ray Bursts due 
their fast rotation. During about 1.5 million 
years, in the very early phases of the globu- 
lar cluster (typically with ages between 3.5 and 
5 million years) one expects a rate of GRB 
of about 2 per 1000 years per globular clus- 
ter, or 2 per 10 years when the whole system 
of globular clusters of a galaxy is surveyed. In 
case the slope of the IMF is 1.35, the expected 
rates are lowered by about three orders of mag- 
nitudes. Typically the GRB rate would be be- 
tween 0.002 and 0.004 per 1000 years and per 
globular clusters. 

4.3. He-rich stars and the UV-upturn in 
ellipticals 

Interestingly, all ellipticals show the UV- 
uptum phenomenon, i.e. an UV excess in 
their spectra. These excesses are caused by 
an old population of hot helium-burning stars. 
Among the models which have been prop osed 
to exp l ain this UV-upturn, let us mention iLed 
(11994 : iPark & Led (1 19971) who propose that 
the blue core-helium-burning stars, originate 
from the low-metallicity tail of a stellar pop- 
ulation with a wide metallici t y dist ribution. 
iBressan et aP (Il994l) : lYi et al.l (1 19971) on the 
contrary, assumes that the UV-upturn is caused 
by metal-rich stars that lose their envelope near 
the tip of t he firs t giant branch. More recently, 
iHan et al propose that these stars lost 

their envelope because of binary interactions. 

As seen above. He-rich stars seem to 
be present in globular clusters (already on 
the Main Sequence). As in elliptical galax- 
ies, some globular clusters also present an 
extended blue horizontal branch. In globular 
clusters this hot branch of helium-burning stars 
is believed to arise from the He-rich stellar 
population. Could it be also the case in ellip- 
tical galaxies? Or said in other words, do ellip- 
ticals also contain a population of non-evolved 
He-rich low mass stars, whose origin may be 
linked with a previous generations of fast ro- 
tating massive stars? 

In the monolithic formation scenario, ellip- 
tical galaxies underwent at their birth a very 
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strong burst of star formation. If there is some 
link between strong star formation episodes 
and the presence of very fast rotators, then one 
would expect in elliptical galaxies a great num- 
ber of massive fast rotators. Would these fast 
rotators be able to produce anomalous stars 
in ellipticals? A big difference between glob- 
ular clusters and ellipticals is the depth of the 
gravitational potential well. In ellipticals, the 
slow and fast winds and the supernovae ejecta 
are retained (at least in part). Thus the cir- 
cumstances are less favorable to the appari- 
tion of stars, made of a mixture of interstel- 
lar medium and of only H-burning product^ 
However this might not be correct and the ar- 
gument is the following: in our scenario for 
globular clusters, we must assume that the star 
formation occurs rapidly in the vicinity of the 
massive stars when the slow winds is released. 
If not, the fast winds and the supernovae ex- 
plosions would remove this interstellar mate- 
rial. Interestingly, clumps or protostars are ob- 
serve d in the disk of the (Herbig) Be star MWC 
1080 (IWang et al.l2007l) . indicating at least that 
low mass star formation can occur at a very 
early stage in an equatorial circumstellar disk. 
Thus in ellipticals, even if fast winds and su- 
pernova ejecta would be retained, provided the 
formation of clumps in equatorial disk occur 
rapidly we still could have stars with strong he- 
lium enrichments and bearing the signature of 
CNO-processed material. Thus a similar sce- 
nario than the one suggested above for ex- 
plaining the anticorrelation in globular clus- 
ters may have occurred in elliptical galaxies. 
This may have interesting consequences on the 
early phases of the photometric and chemical 
evolution of the massive elliptical galaxies. 
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